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12.  AfITRAC  T 


e  results  of  a  computational  study  or  one  aspect  of  conden- 
sation  in  the  exhaust  plume  consisting  of  H2O  and  Hz  of  an  Apollo  rocket 
engine  are  reported.  Saturation  of  water  vapor  occurs  at  about  225^K  in 
the  plume.  Hydrogen  remains  unsaturated  in  the  given  flow  field.  Under 
the  assumption  of  steady  state  nucleation,  onset  of  H2O  condensation  by 
homogeneous  nucleation  is  estimated  to  occur  at  x=47,  40,  and  13  m,  i.e., 
about  45,  38,  and  13  m  downstream  of  the  nozzle  exit  plane,  along  9=0. IT, 
0.3,  and  0.8  stream  tubes,  respectively.  It  is  further  estimated  that 
at  xslSO  m  along  the  center  stream  tube  a  cloud  of  10*  *  H2O  droplets/cm* 
is  formed.  The  mass  mean  radius  of  this  cloud  of  droplets  is  about  17  a, 
ana  there  are  only  about  1/100  as  many  droplets  with  r=24  %  or  21  A  as 
with  r=17  A. 

It  must  be  emphasized  that  the  fundamental  assumption  of  steady  stat 
nucleation,  on  which  the  above  conclusions  are  based,  appears  to  be 
questionable  for  this  problem.  The  flow  is  too  fast  for  the  steady  state 
nucleation  rate  to  be  attained*  Consequently,  H2O  condensation  by  homo¬ 
geneous  nucleation  may  be  delayed  or  even  inhibited  in  this  exhaust  plume 

However,  we  must  remember  that  homogeneous  nucleation  is  by  no  means 
the  only  mechanism  that  can  induce  H2O  condensation.  Earlier  condensatio 
may  be  initiated  by  the  processes  of  heterogeneous  and  binary  nucleation. 
These  processes  are  being  studied  at  present. 
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I  INTRODUCTION 


The  results  of  a  computational  study  of  one  aspect  of 
condensation  phenomena  in  the  high  altitude  exhaust  plume 
of  an  Apollo  rocket  engine  are  presented  here.  The  exhaust 
gas  consists  of  water  vapor  and  hydrogen.  The  objective 
of  this  study  is  to  determine  whether  any  of  the  species 
is  likely  to  condense,  and  if  sc,  at  what  location  in  the 
plume  the  onset  of  condensation  may  be  expected. 

Only  the  possibilities  of  condensation  by  the  process 
of  homogeneous  nucleation  have  been  studied.  The  effects 
of  the  other  processes  which  may  induce  condensation,  namely, 
binary  nucleation  and  heterogeneous  nucleation  have  not  yet 
been  studied  for  this  investigation. 

In  addition  to  nucleation  the  average  growth  of  droplets 
was  calculated  beyond  the  point  of  onset  of  condensation 
and  finally  a  rough  estimate  of  the  droplet  size  distribution 
far  downstream  of  the  onset  point  has  been  made. 

The  general  method  of  calculation  is  the  same  as  that 
employed  in  an  earlier  ~tudy  of  the  exhaust  plume  of  a 
different  rocket  engine  (Wu  1972).*  The  method  will  there¬ 
fore  not  be  reiterated  here  in  detail. 

Tables  I  -  IV  summarize  the  flow  variables  of  the 

exhaust  plume  in  three  stream  tubes  designated  as  d  ~  0.0, 

0.3,  and  0.8  radians  as  transmitted  to  us  by  Teare  (1972) 

*  Note  references  at  the  end  as  listed  alphabetically  by 
senior  author.  Moreover,  because  of  the  constant  reference 
to  the  earlier  report  (Wu  1972),  it  is  reproduced  here,  by 
permission  of  AeroChem  Research  Laboratories,  Inc..  Princeton 
N.J.,  as  a  supplement  to  this  report. 


of  AVCO-Everett  Research  Laboratory.  The  first  data  point 
given  for  each  stream  tube  corresponds  to  the  location  of 
the  exit  of  the  rocket  nozzle.  As  before,  it  was  suggested 
to  us  that  the  flow  field  is  chemically  frozen  in  the  region 
of  the  given  data.  Hence  the  law  of  mass  conservation  can 
be  applied  to  ^0  and  independently.  In  the  absence  of 
appreciable  condensation,  the  partial  pressures  of  the 
species  are  easily  computed  by  taking  the  mole  fractions 
of  the  species  to  be  constant.  Figure  (1)  shows  the  variation 
of  the  total  and  partial  pressures  of  the  species  with 
temperature  along  the  center  stream  tube  (0  =  0.0).  The 
location  on  the  stream  tube  is  not  marked.  The  pressure  - 
temperature  functions  for  the  other  two  stream  tubes  are 
practically  indistinguishable  from  those  in  Fig.  (1).  The 
saturation  vapor  pressure,  p^  Q(T),  water  a3-so 
pitted  in  Fig.  (1).  The  intersection  of  the  saturation 
line  with  the  partial  pressure  curve  marks  the  state  at 
which  the  water  vapor  in  the  exhaust  becomes  saturated. 

This  location  occurs  at  x  =  3220  cm  on  the  center  stream 
tube  which  is  about  3050  cm  beyond  the  exit  plane.  The 
saturation  conditions  on  all  three  stream  tubes  are  summarized 
in  Table  V.  The  lowest  temperature  in  the  given  flow  field 
is  69°K.  Since  the  critical  temperature  of  hydrogen,  43°K, 
is  still  lower, this  vapor  remains  unsaturated  in  the  exhaust. 

The  effective  molecular  weight  and  ratio  of  specific 
heats  of  the  exhaust  gas,  computed  under  the  assumptions 
of  chemically  frozen,  isentrdpic  flow,  are 


U  s  13.102 


and 


1.217  <  y  <  2.351  . 

As  in  the  previous  case,  the  value  of  y  varies  somewhat 
along  the  stream  tube.  However,  for  69  <  T  <  22 5° K,  y 
remains  practically  coni,  cant  at  1.351  for  all  three  stream 
tubes.  Since  the  condensation  calculations  would  involve 
only  the  flow  region  in  which  the  gas  temperature  falls  in 
this  limited  range,  it  was  decided  to  use  the  value  y  ~  1.351 
in  the  calculation  of  the  entire  flow  field  until  the 
point  where  appreciable  condensation  has  noticeably  altered 
the  gas  phase  composition. 

With  the  molecular  weight  and  ratio  of  specific  heats 
known,  the  supply  conditions*  of  a  "fictitious  reservoir", 
p0  and  T0 .  and  the  area  distribution*  of  a  "fictitious 
nozzle"  are  easily  computed  using  the  data  in  Tables  I  -  IV 
and  the  equations  of  one-dimensiona2.  isentropic  flow  of  ideal 
gases.  The  values  of  these  parameters  are  listed  in  Table  VI, 
and  they  are  used  as  the  input  to  the  computer  programs. 

It  should  be  pointed  out  that  the  use  of  these  "fictitious 
parameters"  is  only  practical  for  calculations.  The  parameters 
associated  with  the  fictitious  nozzle  and  reservoir  used 
here  have  no  bearing  on  the  actual  rocket  nozzle  and  supply 
chamber  conditions.  It  is  noted  that  a  simple  square  lav; 

A/A*  =  C  x2  (1) 

*  These  quantities  are  calculated  to  facilitate  subsequent 
computations,  see  pp.  9-11  of  Wu  (1972)  for  details. 
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is  used  ho  describe  the  area  variation  along  each  stream 
tube.  This  is  an  accurate  correlation  for  the  part  of 
the  stream  tube  where  the  condensation  calculations  will 
be  made,  as  seen  in  Figure  (2). 

The  mass  fraction  of  the  condensate  is  computed  by 
integrating  the  equation  (see,  e.g.  Wegener  and  Mack  1958) 

4  TT  p  X 

g(x)  =  —5“  f  J(v)A(?)r3(x,?)d?  .  (2) 


Here,  r(x,£)  is  the  final  radius  of  a  droplet  initially  formed 
a.t  C  by  homogeneous  nucleation  and  having  undergone  growth 
before  reaching  the  poirt  x  downstream  of  5.  This  radial 
growth  is  expressed  formally  by 


r(x,?)  =  r*CO  + 


dr 

Hx 


dC  , 


(3) 


where  dr/dx  is  the  ..roplet  growth  rate.  In  addition,  J(c) 
is  the  homogeneous  nucleation  rate  of  critical  clusters  taken 
to  be 


J 


Jclass  ’ 


(4) 


where  T  is  an  empirical  correction  factor  and  the  value  of 
Jciass  "t^ie  rate  to  be  computed  from  the  so-called  classical 
theory  of  homogeneous  nucleation.  With  a  proper  choice  of  T, 
we  avoid  the  difficulties  involved  in  calculating  J,  caused 
by  our  lack  of  knowledge  of  the  surface  tension  of  small 


systems  and  the  uncertainties  in  nucleation  theory.  The 
droplet  growth  rate  is  computed  from  the  free  molecular 
equation  as  discussed  in  the  previous  report  (Wu  1972), 
since  again  the  droplet  radii  are  much  smaller  than  the 
vapor  mean  free  path. 


II  RESULTS  AND  DISCUSSIONS 


As  a  direct  check  on  the  input  data  to  the  computer 
programs,  the  computed  pressure  -  temperature  relation  along 
the  center  stream  tube  near  the  saturation  point  of  water 
vapor,  together  with  the  initially  given  data  points  are 
plotted  in  Figure  3.  The  agreement  is  seen  to  be  good. 

Figure  (4)  shows  the  accummulation  of  H20  condensate 
mass  fraction  g  =  mc/m  along  the  center  stream  tube  as  a 
function  of  the  empirical  correction  factor  F  here  indicated 
in  powers  of  ten.  The  qualitative  feature  of  an  onset  of 
condensation  is  clearly  exhibited  in  these  curves  with  the 
appearance  of  the  condensate  as  a  sudden  process .  This 
qualitative  feature  is  well  verified  in  many  observations 
and  it  does  not  depend  on  the  choice  of  F.  The  flow  time 
of  fluid  elements  originating  at  H20  saturation  is  shown 
by  the  time  of  flight  scale  at  the  bottom  of  this  figure. 

Defining  the  onset  of  condensation*  as  the  point  where  the 

-3 

mass  fraction  of  condensate  g  =  10  ,  we  plot  the  calculated 

point  of  onset  in  the  center  stream  tube  in  Fig.  (5)  again 
as  a  function  of  the  correction  factor  T.  Similar  results 
for  the  flow  in  the  other  two  stream  rubes  are  shown  in 
Figures  (6)  and  (7). 

The  lower  limit  of  T  is  likely  to  be  of  the  order  one, 
which  corresponds  to  the  homogeneous  nucleation  rate  given 
by  the  classical  theory. 

*  The  choice  of  this  definition  1ms  been  discussed  in 
detail  in  the  pr  vious  report  (Wu  1972). 


On  the  other  limit,  the  theoretical  value  of  T  is  estimated 
to  be 
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^tr  ^rot  _  IQ 14 


*  .  3 

where  q  ,  the  replacement  factor,is  given  the  value  10 

X  v  U 

as  proposed  by  Lothe  and  Pound  (1962)  in  their  statistical 

mechanical  model.  The  translational  and  rotational  partition 
.  ~  *  * 

functions,  q  and  q  t,  are  computed  from  well  known  formulas 
from  statistical  mechanics  (e.g.  Wegener  and  Parlange  1970). 
This  latter  value  of  P  is  much  higher  than  the  empirical 
estimates  discrepancy  between  Lothe  and  Pound’s  theory  and 
experiments  found  before  at  higher  temperatures .  Other 

ft 

models  for  q  (e.g.  Dunning  1965)  yield  values  of  T  that 

j. 

are  much  lower.  However,  none  of  the  theories  proposed  are 

yet  in  a  state  to  be  applicable  to  engineering  estimates. 

14 

We  note  in  Figure  (5)  that  the  range  1  <  r -theory  <  10 
corresponds  to  the  range  41  <  xonset  <  55  m  calculated  for 
steady-state  homogeneous  nucleation  theories  in  the  center 
stream  tube.  In  view  of  this  uncertainty,  we  shall  use 
empirical  knowledge  gained  in  previous  experiments  with 
condensation  in  nozzles. 

For  this  work, therefore, the  choice  T  -  10°,  a  -  a.  = 

ice 

96  dyne/cm  and  a  =  1.0  has  again  been  made  to  estimate  the 
onset  of  condensation  of  HjO.  The  motivation  for  this 
assumption  for  the  empirical  correction  has  been  further 
elaborated  in  Wu  (1972).  The  choice  of  these  values  was 
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based  on  experiments  where  the  mole  fraction  of  the  condensing 
vapor  y  <<  1.  In  this  nozzle,  however,  we  have  q  -  0.6929, 
a  high  value  that  necessitates  an  examination  of  the  non- 
isothermal  kinetics  of  nucleation. 

The  so-called  ncn-isothermal  effect  in  homogeneous 
nucleation  was  first  noted  by  Kantrowitz  (1951),  treated 
in  detail  by  Feder  et  al.  (1966),  and  applied  to  condensation 
in  steam  nozzle  flows  by  Barschdorff  et  al.  (1972).  It  deals 
with  the  fact  that  the  clusters  are  generally  not  at  the 
same  temperature  as  the  surrounding  gas ,  and  that  the  latent 
heat  of  condensation  added  to  a  cluster  when  it  is  growing 
towards  the  critical  size  .-»>ust  be  dissipated  to  the 
surrounding  gaseous  medium.  This  heat  transfer  process  is 
coupled  with  the  growth  of  the  clusters  and  results  in 
a  slower  growth  rate  of  clusters  to  their  critical  size. 

The  nucleation  rate  will  in  turn  be  retarded  due  to  the 
transfer  of  latent  heat.  Hence,  one  expects  the  non-isothermal 
nucleation  rate  to  be  lower  than  the  isothermal  nucleation 
rate  given  by  the  classical  theory.  When  the  vapor  is 
carried  in  a  large  quantity  of  inert  diluent  which  participates 
in  heat  transfer,  the  non-isothermal  nucleation  rate  may 
approach  the  isothermal  rate.  The  non-isothermal  nucleation 
theory  of  Feder  et  al.  (1966)  has  recently  been  examined 
by  Wilemski  (1973)  who  found  that  the  assumptions  made  by 
Feder  et  al.  are  indeed  reasonable  and  their  theory  is 
applicable  to  this  work.  A  simple  calculation  using  the 
results  of  Feder  et  al.  (1966)  reveals  that 
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noniso 


b  +q‘ 


^iso  ’ 


we  have 


°*02  <  ~2 — 7  <  0.024  , 

b  +q 


for  3250  <  x  <  4680  cm  i„  the  center  stream  tube.  The 
quantities  b  and  q  are  defined  in  terms  of  the  thermodynamic 
state  of  the  supersaturated  vapor  and  the  inert  carrier 
gas ,  and  applicable  expressions  for  these  quantities  can 
be  found  in  Wu  (1072).  Hence,  the  uon-isothermal  nucleation 
factor  is  found  to  be  nearly  a  constant  of  the  order  of  10~2. 
In  other  words .  the  non-isothermal  nucleation  is  only  about  ’ 
1/100  of  the  standard  one.  As  before,  we  estimate  that 
the  proper  value  of  r  for  the  problem  being  considered 
should  lie  between  lo"  and  106.  Moreover,  it  is  implied 
throughout  that  the  assumptions  underlying  the  steady-state 
theory  of  homogeneous  nucleation  are  valid.  From  these 
calculations,  we  finally  find  that  condensation  of  water 
vapor  may  be  expected  in  the  plume  at 


X  5  47  m  for  6  :  0.0  stream  tube, 

x  =  40  m  for  6  =  0.3  stream  tube,  and 

x  s  13  m  for  6  =  0.8  stream  tube. 

These  x-coordinates  correspond  roughly  to  45,  38,  and  13  m 

downstream  of  the  nozzle  exit  plane  along  9  ,  0.0,  0.3, 


and  0.8  stream  tubes  ^respectively. 


The  thermodynamic  conditions  at  the  onset  of  conden- 

5 

sation  computed  for  the  three  stream  tubes  assuming  F  =  10 
are  summarized  in  Table  VII.  Moreover,  the  onset  condition 
in  the  center  stream  tube  is  compared  with  available  experi¬ 
mental  results  in  Fig.  (8).  It  is  again  noted  that  the 
temperature  and  vapor  pressure  at  the  onset  of  condensation 
are  still  much  lower  than  the  range  of  states  for  which 
experiments  are  available.  However,  the  pressures  are 
higher  than  those  found  in  the  previous  case  (Wu  1972'  where 

_3 

Ph0»10  torr  at  the  onset  with  T  -  150°K.  It  is  our 
opinion,  however,  that  extrapolation  of  our  empirical 
knowledge  on  F  to  the  temperatures  and  pressures  occurring 
in  this  problem  is  not  unreasonable.  We  shall  see  later 
that  there  are  sources  for  more  serious  uncertainties ,  which 
may  outweigh  the  uncertainties  stated  so  far. 

Under  the  assumption  that  the  stream  tube  area  variation 
is  not  affected  by  the  latent  heat  addition  to  the  gas 
stream  resulting  from  condensation,  calculations  of  droplet 
growth  beyond  the  onset  point  can  be  made  with  the  available 
computer  programs .  The  results  of  such  a  calculation  are 
shown  in  Figure  9.  Here,  r  is  the  mass  mean  radius  of  all 
the  droplets  present  ,  i.e.,  the  spherical  radius 

resulting  from  evenly  distributing  the  total  mass  of 
condensate  over  all  the  droplets  present.  Thus, 


9 


1/3 
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where  p  is  the  density  of  the  vapor-carrier-  gas  mixture  and 
N  is  the  total  number  of  drop.lets  present.  The  droplet 
size  distribution  was  calculated  by  following  the  growth 
process  of  particular  classes  of  drople  :s  as  they  are  carried 
downstream  by  the  gas  flow.  Each  class  consists  of  droplets 
formed  at  a  particular  location,  C,  along  the  stream  tube. 
Hence,  an  observer  located  somewhere,  x,  downstream  of  this 
location  will  see  a  total  of  J(?)  A{?)  d£  droplets  in  this 
class  flying  by  per  unit  time.  All  of  the  droplets  in 
this  class  will  have  radii  between  r  and  r  +  dr,  where  r 
is  given  by  Eq.  (3),  and  dr  is  related  to  d£*  Because  the 
nucleation  zone  is  quite  narrow,  the  cross-sectional  area 
of  the  stream  tube  does  not  change  much  in  this  particular 
zone  and  the  number  of  droplets  in  different  classes  will 
be  roughly  proportional  to  the  nucleation  rates,  J(c),  at 
the  respective  "birth  places"  of  the  classes  of  droplets. 

Hence,  the  class  having  the  largest  number  of  droplets 
must  be  formed  at  the  location  where  J  is  at  its  maximum. 

In  the  center  stream  line,  with  the  assumption  r  =  10^, 

occurs  at  x  =  47C0  cm.  The  radii  of  droplets  formed 

max 

at  x  =  4700  cm  [i.e.  r(4700,  x)]  practically  coincides 
with  r  forx  >  4900  cm.  This  means  that  the  "peak"  in  the 
droplet  size  distribution  coincides  with  the  mass  mean  radius 
of  all  droplets  present.  This  qualitative  feature  is  well 
established  by  laser  light  scattering  experiments  with  water 
vapor  (Wegener  and  Stein  1969)  and  with  ethyl  alcohol  by 
Clumpner  (1971). 

Four  more  curves,  r(4440,  x) ,  r(4C40,  x) ,  r(4880,  x) 
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and  i'(4940,  x)  are  also  shown  in  Fig.  9.  The  nucleation 
rates  calculated  at  these  points  satisfy  the  following 
relation: 

J (4700)  :  JC4880)  :  J(4940)  =  J(4700)  :  JC4540)  :  JC444Q) 

=  1  :  0.1  :  0.01  . 


Consequently,  one  estimates  that  if  the  assumptions  under¬ 
lying  these  calculations  are  valid,  there  will  be  ten  times 
fewer  droplets  in  the  classes  r(4880,  x)  and  r(4540,  x> 
than  in  the  class  r(4700,  x)  which  has  the  largest  number 
of  droplets.  Similarly,  there  will  be  one  hundred  times 
few«r  droplets  in  the  classes  r(4940,  x)  and  r(4440,  x)  then 
the  class  r(4700,  x).  In  other  words,  we  estimate  that  at 
x  s  150  n  in  the  center  stream  tube,  the  mass  mean  radius 
of  all  the  droplets  is  17  X  and  this  happens  to  be  the 
"most  probable"  size  of  the  droplets.  An  H20  droplet  of  this 
size  contains  about  600  molecules.  At  this  size,  macroscopic 
notions  like  surface  tension  are  probably  meaningful. 
Furthermore,  we  estimate  that  the  numbers  of  droplets  with 
r  =  15  X  and  20  X  are  only  1/10  of  that  of  17  X  and  those 
of  14  A  and  21  X  are  only  1/100  of  the  number  of  droplets 
with  r  =  17  X.  The  total  number  density  of  droplets  of  all 
sizes,  computed  from 


> 


N(x) 


J(C)A(OdC 
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10  -3 

is  found  to  be  in  the  order  of  1C  cm  at  x  =  150  m. 

We  must  now  discuss  the  weakest  link  in  the  above 
analysis;  all  the  above  calculations  are  based  on  the 
fundamental  assumption  of  a  steady-state  homogeneous 
nucleaticn  process.  A  careful  examination  of  the  thermodynamic 
state  and  flow  variables  in  the  supersaturated  region  reveals  , 
however,  that  this  assumption  may  not  be  valid.  The  flow 

5 

time,  t^.,  from  saturation  to  onset  calculated  for  T  =  10  ‘ 

is  about  3  msec  along  the  center  stream  tube  as  the  bottom 

scale  in  Fig.  4  indicates.  The  build-up  time,  t  ,  for 

steady-state  homogeneous  nucleation  under  the  existing 

thermodynamic  conditions,  calculated  using  the  ideas  of 

Kantrowitz  (1951),  is  about  14  msec.  The  appropriate 

Damkdhler  parameter,  D  =  T^/Tp  is  therefore  D  a  0.2  or 

D  <  1  in  the  exhaust  plume  being  considered.  A  typical 

value  for  D  found  for  condensation  of  ethanol  in  the  nozzle 

2 

used  by  Wegener  et  al.  (1972)  is  of  the  order  10  .  The 
condition  D  >  1  must  be  satisfied  in  order  that  the  steady 
state  cluster  distribution  may  be  established  when  the 
thermodynamic  state  of  the  vapor  is  changing.  As  in  the 
previous  case,  we  found  here  that  there  is  probably  not 
enough  time  for  the  steady  state  nucleation  rate  to  be 
attained.  The  actual  homogeneous  nucleation  rate  may 
therefore  be  lower  than  the  steady-state  rate  used  in  the 
calculations  here.  Hence,  a  calculation  involving  non¬ 
steady  homogeneous  nucleation  "ate  may  yield  onset  of 
condensation  further  downstream  from  that  found  in  this 
work,  or  it  may  predict  no  condensation  at  all.  At  any 
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p- 

rate,  condensation  by  homogeneous  nucleation  would  most 
[  likely  not  occur  up- -‘•ream  of  the  onset  locations  found  here. 
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III  SUMMARY  AND  CONCLUSIONS 


Possible  condensation  by  homogeneous  nucleation  is 
investigated  for  the  exhaust  plume  of  an  Apollo  rocket 
engine.  The  exhaust  studied  consists  of  a  chemically  frozen 
mixture  of  HjO  and  H2.  Water  vapor  becomes  saturated  at 
aiout  225°X  in  the  exhaust  and  subsequent  expansions  to 
lower  temperatures  may  lead  to  its  condensation.  Hydrogen 
remains  unsaturated  in  the  region  for  which  the  flow  field 
is  given. 

Onset  of  water  vapor  condensation  induced  by  homogeneous 

nucleation  has  been  calculated  assuming  steady-state 

nucleation  theory  for  non-isothermal  cluster  growth  and  using 

empirical  information  gathered  in  numerous  previous  studies 

at  higher  temperatures  and  pressures.  The  results  of  this 

calculation  show  that  H20  condensation  will  take  place  at  about 

170°K,  which  is  approximately  45,  38,  and  13  m  downstream 

of  the  nozzle  exit  plane  for  P  =  0.0,  0.3,  and  0.8  radian 

stream  tubes,  respectively.  In  the  center  stream  tube 

(0  =  0.0),  calculations  of  the  condensation  process  beyond 

the  onset  point,  using  the  stream  tube  area  distribution 

of  nonccndensing  flow  lead  to  an  estimate  of  the  droplet 

size  distribution.  It  is  estimated  that  at  x  =  150  m  in 

10  0 

the  center  stream  tube,  about  10  droplets/cm"  are  formed 
by  homogeneous  nucleation.  The  mass  mean  radius  of  this 
cloud  of  droplets  is  about  17  X.  This  is  also  the  size  of 
the  largest  number  of  droplets,  i.e.  the  size  distribution 
is  quite  "narrow1*.  There  are  only  1/10  as  many  droplets 
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with  r  =  IS  X  and  20  X  as  droplets  with  r  =  17  X,  and  1/1G0 
with  r  =  14  X  and  21  X.  There  are  about  600  molecules 
in  a  droplet  with  r  -  17  X. 

It  must  be  emphasized  that  the  above  conclusions  are 
derived  from  the  crucial  assumption  of  steady-state 
homogeneous  nucleation.  A  further  examination  of  the  flow 
field  reveals  that  this  assumption  does  not  appear  to  be 
valid  for  the  problem  at  hand.  The  build-up  time  for  steady 
state  homogeneous  nucleation  is  estimated  to  be  14  msec 
while  the  relevant  flow  time  is  only  about  3  msec .  Thus , 
it  is  seen  that  the  flow  is  vrobably  too  rapid  for  the 
steady  state  nucleation  rate  to  be  achieved.  Consequently, 

1^0  condensation  induced  by  homogeneous  nucleation,  as 
summarized  in  the  preceeding  paragraph  may  occur  at  locations 
farther  downstream  than  estimated  or  it  may  hot  cccur  at  all. 

On  the  other  hand,  condensation  by  homogeneous  nucleation 
cannot  be  ruled  out  at  the  present  time  since  the  steady 
sta-'  i  theory  provides  only  one  limit  for  its  occurrence. 

We  must  finally  emphasize  that  homogeneous  nucleation  is  not 
the  only  mechanism  capable  of  initiating  HjO  condensation. 

The  processes  of  heterogeneous  and  binary  nucieation  are 
also  possible  ar.d  these  mechanisms  both  faver  early  phase 
change.  Work  in  these  directions  is  currently  being  carried  out 
at  Yale  since  it  is  required  to  assess  the  probability 
of  condensation  caused  by  heterogeneous  or  binary  nucleation. 
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LIST  or  SYMBOLS 


a  speed  of  sound. 

A  cross  sectional  area, 

g  mass  fraction  of  condensate. 

J  nucleation  rate. 

«s  mass  flow  rate  in  the  nozzle. 

M  Mach  number . 

n  number  of  molecules  in  a  cluster. 

N  number  density  of  clusters, 
p  pressure. 

p—  saturation  vapor  pressure  ever  a  flat  surface. 

q  replacement  factor, 

rep 

5  ,  q  translational  and  rotational  partition  functions - 

tr  rot 

r  radius . 

s  saturation  ratio  (  =  p^/p^) • 

t  tine . 

T  temperature . 

u  velocity. 

x  coordinate  along  flow  direction, 

y  nole  fraction, 

a  n=ss  ^ccossodaticn  coefficient. 

Y  ratio  of  specific  heats, 

f  empirical  factor, 

p  molecular  weight, 

p  mass  density. 

o  surface  tension  (or  surface  free  energy). 
t p  characteristic  flow  tine. 
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Tr  relaxation  time  for  steady-state  nucleation. 
w  mass  fraction  of  vapor. 

Subscripts 
c  condensate, 

class  classical, 

i  inert  carrier  gas. 

noniso  nonisothermal . 

o  initial  supply  (stagnation)  conditions, 

r  a  cluster  of  radius  r. 

v  condensable  vapor. 

00  flat  surface. 


Superscript 

cluster  of  critical  size,  or  nozzle  throat. 
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Table  I 

SPECIFICATIONS  FOR  CONDENSATION  CALCULATION  II 

Initial  conditions  for  3  stream  tubes: 

T  =  1417°K 

p  =  125.4  torr 

p  =  1.860  x  10"5  g/cm3 

u  =  4.170  x  10-3  cm/sec 

2 

Area  =  3.010  cm 

[HjO]  -  *6929  -j  moie  fractions 
[H23  =  .3071  J 

All  three  stream  tubes  have  same  initial  gas  properties  at 

2  4  2 

area  A  =  3.01  cm  but  expansion  to  A  *  3  x  10  cm  occurs 

over  different  axial  distances: 

0=0  169.5  <  x  <  17320  cm 

0=0.3  143.8  i  x  i  14590  cm 

0=0.8  46.79  i  x  t  4793  cm 


Table  II 


x  (cm) 

169.5 

199.7 
235.9 

277.7 

325.7 
396.2 

480.5 

603.5 
756.0 
1017 
1417 
2117 
3788 
10100 
17320 


8  =  0.0  Radians 


T 

P 

u 

A 

125.4 

1417 

1 . 86xlG~5 

4.17xl05 

3.010 

83.37 

1311 

1.33 

4.254 

4.108 

54.51 

1202 

9 . 53x10*"^ 

4.332 

5.655 

35.83 

1101 

6.849 

4.401 

7.745 

23.78 

1003 

4.9*8 

4.463 

1.055x10 

14.37 

900.5 

3.353 

4.531 

15.37 

8.717 

803.3 

2.281 

4. 590 

22.30 

4.805 

698.4 

1.446 

4.652 

34.72 

2.655 

605.4 

9.214xl0"7 

4.704 

53.86 

1.207 

498.5 

5.087 

4.762 

95.38 

0  c  4980 

398.6 

2.424 

4.815 

1.848x10 

0.1690 

302.0 

1.175 

4.864 

4.084 

0.03508 

200.8 

3.67xl0”8 

4.914 

1.294x10 

2.475xl0“3 

100.8 

5.158xl0"8 

4.964 

9.119 

5.765xl0-4 

69.0 

1.755 

4.979 

2.671x10 

X  (cm) 

_E~ 

143. 8 

1*5.4 

169.1 

83.52 

199.5 

54.89 

234.4 

36.15 

274.7 

24.12 

333.9 

14.61 

404.7 

8.877 

508.2 

4.89  5 

o36 . 5 

2.706 

856.6 

1.230 

1193 

0.5078 

1782 

0.1723 

3189 

0.03578 

8307 

2.524xl0"3 

14590 

5.869xl0“4 

T 

P 

1417 

1.86xl0“5 

1310 

1.34‘ 

1202 

9.590xl0“6 

1102 

6.891 

1009 

5.021 

902.7 

3.4PO 

805.5 

2.315 

700.6 

1.468 

607.6 

9.3S9xl0‘7 

500.3 

5.168 

400.2 

2.666 

303.1 

1.194 

201.6 

3.729xl0~8 

101.2 

5.24xl0“9 

69.3 

1.78 

u 

A 

4.17X105 

3.010 

4.253 

4.098 

4.330 

5.622 

4.399 

7.701 

4.461 

1.042X101 

4.529 

1.516 

4.588 

2.198 

4.649 

3,420 

4.702 

5.305 

4.760 

9.490 

4.813 

1.819xl02 

4.862 

4.021 

4.913 

1.274xl03 

4.953 

8.978 

4.973 

2.634xl04 

T 


P 


u 


A 


x  (cm) 

T 

P 

u 

46.79 

125.4 

1417 

1.86xl0“5 

4.17x10 

54.78 

84.51 

1309 

1.356 

4.251 

64.17 

56.73 

1206 

9.885xl0“6 

4.325 

75.06 

38.12 

1108 

7.226 

4.392 

87.68 

25.64 

1017 

5.296 

4.453 

10G.3 

15.64 

911.4 

3.605 

4.520 

128.6 

9.546 

814.4 

2.463 

4.579 

161.2 

5.287 

709.3 

1.566 

4.641 

201.8 

2.929 

615.6 

9.996xl0”7 

4.694 

261.7 

1.471 

519.9 

5.947 

4.747 

364.4 

0.6075 

416.3 

3.066 

4.802 

544.7 

0.2061 

315.6 

1.372 

4.853 

940.1 

0.04723 

215.4 

4.607xl0"8 

4.904 

2249 

4 .  469x10*" 3 

116.7 

8.050xl0~9 

4.953 

4793 

S.784xl0"4 

68.6 

1.772 

4.976 

3.010 

4.049 

5.460 

7.356 

9.899 

1.433X101 

2.070 

3.212 

4.975 

8.271 

1.5 86x10 2 

3 . 505 

1.033xl03 

5. 856xi03 

2.648xl04 


Table  V 


x(cm) 

T(°K) 
p  (torr) 

Ph2o 

M 

u  (cm/sec) 
A/A* 


Saturation  Conditions 


e  =  o.o 

9  =  0.3 

3225 

2720 

225 

225 

0.054 

0.055 

0.038 

0.038 

11.16 

11.15 

490100 

489700 

9  =  0.8 
860 
229 
0.060 
0.042 

11.14 

493600 


1854 


1845 


1837 
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Table  VI 


Summary  of  Input  Conditions 


0  =  0.0 

CO 

• 

o 

II 

a> 

CO 

• 

o 

II 

<r> 

p0  (torr) 

9249 

9318 

9485 

T0  (°K) 

5146 

5142 

5138 

“»o 

0.9527 

0.9527 

0.9527 

yi 

2.016 

2.016 

2.016 

18.016 

18.016 

18.016 

1.403 

1.403 

1.403 

Yv 

1.332 

1.332 

1.332 

A*  (cm2) 

0.5058 

0.5021 

0.4926 

C  (cm-2) (a) 

1.783xl0-4 

2.495xl0"4 

2.373x10 
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Table  VII 

Computed  conditions  at  onset,  using  r  =  105,  classical 
steady-state  homogeneous  nucleation  theory  with  a  -  a 

00 


O 

• 

o 

II 

CD 

9  =  0.3 

9  =  0.8 

x  (m) 

47.0 

39.6 

13.0 

p  (torr) 

0.0195 

0.0199 

0.0195 

J>„20  (torr) 

0.0135 

0.0138 

0.0135 

T  (°K) 

173 

173 

171 

A/A* 

3940 

3910 

4010 

M 

12.8 

12.8 

12.9 

scr 

"1500 

"1450 

"2000 

r*(X) 

3.6 

3.6 

3.5 

n* 

6 

6 

6 

N  (cm-3) 

"1010 

"1010 

"1010 

. :  i 


500 


IQOi 


x  (cm) 


Figure  2 


the  -  Pound  theory 


Figure  5 


bothe  -  Pound  theory 


x(  cm) 


Figure  6 


Lothe  -  Pound  theory 
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